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I. INTRODUCTION 



Nonzero ground state entropy, So ^ 0, is an important subject in statistical mechanics. One physical 
example is provided by ice, for which So = 0.82 ± 0.05 cal/(K-mole), i.e., Sq/R = 0.41 ± 0.03, where 
R = Nji vog kg p]-0]. A particularly simple model exhibiting ground state entropy without the complication 
of frustration is the g-state Potts antiferromagnet (AF) on a graph G j?J (which may or may not be a 
regular lattice) for q > x(G), where x(A) denotes the chromatic number of G, i.e., the minimum number of 
colors necessary to color the vertices of the lattice such that no two adjacent vertices have the same color. As 
is evident, there is a deep connection with graph theory here, since the zero-temperature partition function 
of the above-mentioned g-state Potts antiferromagnet on a graph G satisfies 

Z(G,q,T = 0) PAF = P(G iq ) (1.1) 

where P(G, q) is the chromatic polynomial expressing the number of ways of coloring the vertices of the 
graph G with q colors such that no two adjacent vertices (connected by a bond of the graph) have the same 
color (|]- . Since the strict mathematical definition of a graph G requires the number of vertices n — v(G) 
to be finite, we shall denote 

lim G = {G} (1.2) 
In this limit, the ground state entropy per vertex (site) is given by 

S = k B lnW({G},q) (1.3) 
where M^({G}, q), the ground state degeneracy per vertex, is 

W({G},q) = lim P{G,q) 1/n (1.4) 

n — >oc 

As noted, in the limit n — > oo, nonzero ground state entropy, Sq({G}, q) > 0, or equivalently, ground state 
degeneracy W({G} 7 q) > 1 generically occurs in the g-state Potts antiferromagnet for sufficiently large q on 
a given lattice G; equivalently, the number of ways of coloring the graph subject to the constraint that no 
two adjacent vertices have the same color grows exponentially rapidly with the number of vertices on the 
graph. 

Given the fact that P(G, q) is a polynomial, there is a natural generalization, which we assume here, of 
the variable q from integer to complex values. IF({G}, q) is an analytic function in the q plane except along 
a continuous locus of points (and at possible isolated points which will not be relevant here). Following the 



terminology in our earlier papers |12|- p8|, we denote this continuous locus as B. In the limit as n — > oo, 
the locus B forms by means of a coalescence of a subset of the zeros of P(G, q) (called chromatic zeros of G) 
p9[ . Because P(G, q) has real (indeed, integer) coefficients, B has the basic property of remaining invariant 
under the replacement q — > q*: 

B{q) = B{q*) (1.5) 

In a series of papers we have calculated and analyzed the loci B for a variety of families of graphs 

and investigated the connections of these loci with the behavior of W({G},q) for the physical values 
q G Z+, which we have also studied by means of rigorous upper and lower bounds and large-g series [ p^Jl^l 
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and Monte Carlo measurements ||12j,|l3|,|l6[ . In many cases the loci B enclose and form the boundaries of 
two or more regions in the complex q plane. Since this is the case for the families of graphs studied here, 
we shall use the abbreviated nomenclature U W boundary" for B. We note that in other families of graphs, 
such as the infinitely long, finite-width homogeneous strip graphs of various lattices (with open boundary 
conditions) studied in Ref. the loci B form arcs that do not enclose any region. The study of these loci 
for various families of graphs also makes a very interesting connection with algebraic geometry. 

With regard to the properties of the ground state entropy of the Potts antiferromagnet on a given (n — ► oo 
limit of a graph) {G}, two of the most important properties of W({G},q) and the associated nonanalytic 
locus B include (i) the maximal finite value of q at which B crosses the real axis in the q plane, which 
we denote q c , as in our previous work; and (ii) the question of whether B(q) is bounded or unbounded 
(noncompact). Property (i) is important because in cases where B encloses regions, the behavior of the 
ground state degeneracy iy({G}, q) changes qualitatively (nonanalytically) as q (taken to be real) decreases 
through q c . p4j . Concerning the importance of property (ii), we first recall that since an obvious upper 
bound on P(G,q) describing the coloring of an n- vertex graph with q colors is P{G,q) < q n , and hence 
M^({G}, q) < q, it is natural to define a reduced function that is bounded as q — > oo: 

W red X{G},q) = q- 1 W{{G},q) (1.6) 



(This function was denoted W r ({G},q) in our Refs. [[13 15|.) If and only if W re d. ({G}, q) is analytic at 
z = 1/q = 0, there exists a Taylor series expansion of this function around this point. Such large-g series 
expansions are a useful means for obtaining approximate values of the ground state degeneracy of the q-state 
Potts antiferromagnet for finite q on regular lattices |25|- |2^ , [l^ , [l7| . In turn, W re d.({G}, q) is nonanalytic 
at 1/q = if and only if the locus B is unbounded in the q plane, i.e., passes through the origin of the 
1/q plane. Thus, families of graphs with unbounded B(q) do not have large-g Taylor series expansions for 



W re d.{{G}, q). In Ref. [|13| we drew attention to this for bipyramid graphs ]23| and in Ref. |15 we presented 
a general method for constructing families of graphs with unbounded, noncompact loci B. We also presented 
a general (necessary and sufficient) algebraic condition such that B is unbounded. Clearly it is important to 
understand better the differences between the families of graphs that yield W re d.({G},q) functions analytic 
at 1/q = and those that produce W re d.({G},q) functions that are nonanalytic at 1/q — 0. Since B forms 
by merging of chromatic zeros of G as v(G) — > oo, a necessary condition for B to be noncompact in the q 
plane, extending infinitely far from the origin, is that as v(G) — ► oo, the magnitudes \q\ of some chromatic 
zeros must increase without bound. That this is not a sufficient condition is illustrated by the chromatic 
zeros of the complete graph K p psj : these occur at q = 0,1,.., p— 1 so that, as p — > oo, the magnitude of 
the largest chromatic zero grows without bound, but the locus B for this family is not only not noncompact, 
it is trivial, consisting of the empty set. 

In the present paper, continuing our earlier studies |T^,|l5|], we shall address this problem. We shall 
present several additional methods for constructing families of graphs with unbounded B. We shall use these 
to obtain exact calculations of the corresponding chromatic polynomials, W functions, and boundaries B. In 
particular, we shall generate an infinite number of such families via the method of homeomorphic expansion. 
Homeomorphic classes of graphs have been of continuing interest in graph theory [ pT[ ; for some recent 
theorems in a somewhat different direction from the present work, see Ref. [^9| . The study of these families 
is valuable because it enables one gain further insight into the origin of the nonanalyticity of W r ({G}, q) at 
1/q = and the consequent nonexistence of large-g series expansions for the ground state degeneracy of the 
Potts antiferromagnet on such graphs. Here we shall give detailed analyses of families with loci B that have 
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a simple structure in the vicinity of 1/q = 0. In a companion paper, we shall present results for families of 
graphs with loci B that are more complicated in the vicinity of 1/q = 0. 

Before proceeding, we recall two subtleties in the definition of W({G},q) when q is not a positive integer 
p3| . First, for certain ranges of real q, P(G 1 q) can be negative, and, of course, when q is complex, so is 
P(G, q) in general; in these cases it is not obvious, a priori, which of the n roots 

P{G,q) 1/n = {\p(G,q)\ 1/n e ims{p(G ' q))/n e 27Tir,n } , r = 0, 1, ...,n - 1 (1.7) 



to choose in eq. (1A). Consider the function VF({G}, q) defined via eq. ( |l.4| ) starting with q on the positive 



real axis where P(G, q) > 0, and consider the maximal region in the complex q plane that can be reached 



by analytic continuation of this function. We denote this region as R\. Clearly, the phase choice in (|l.7| ) 



for q 6 Ri is that given by r = 0. For families of graphs {G} where there are regions Ri of analyticity of 



W({G},q) that are not analytically connected with Ri, there is no canonical choice of phase in eq. (1.7) 
and hence it is only possible to determine the magnitude ^({G}, q)\ unambiguously. 

A second subtlety in the definition of W({G},q) concerns the fact that at certain special points q s , the 



following limits do not commute ||13| (for any choice of r in eq. (1.7)): 

lim lim P{G,q) 1/n ^ lim lim P(G,q) 1/n (1.8) 
The set q s includes points where P(G, q) has zeros whose multiplicity does not scale like the number of 



vertices |13| (typically these are simple zeros). Following our earlier work 13 1, we define 



W({G},q s ) = W({G},q s ) Dqn = lim lim P(G,r 7 ) 1 /" (1.9) 

q — >q s n — >OG 

This definition has the advantage of maintaining the analyticity of W({G},q) at the special points q s . 
A general form for the chromatic polynomial of an n-vertex graph G is 

N a 

P(G n ,q) = co^ + gc^K-te)^ (1.10) 

where 

t(n)=hn + t (1.11) 

and Cj(q) and aj(q) are certain functions of q. Here the cij(q) and Cj^ (q) are independent of n, while c (g) 
may contain n-dependent terms, such as (—1)", but does not grow with n like (const.) 71 . Obviously, the 
reality of P(G, q) for real q implies that Cj(q) and cij(q) are real for real q. The condition that B does not 
extend to infinite distance from the origin in the q plane is equivalent to the condition that for sufficiently 



large |g|, there is one leading term aj(q) in eq. (1.10). Here we recall that "leading term a^(g) at a point 
q" was defined in Ref. |13| as a term satisfying |a£(g)| > 1 and |a£(g)| > |aj(g)| for j ^ i. (If the cq term is 
absent and N a = 1, then the sole a\(q) may be considered to be leading even if |aj(g)| < 1.) In the limit 
as n — > oo (i.e., for fixed p, r — > oo), the leading term ai in a given region determines the limiting function 
W, with \W\ = \ai\ l1 and the boundary B occurs where there is a nonanalytic change in W as it switches 



between being determined by different leading terms ae in eq. (1.1C). 

Since for the families of graphs studied here, the boundary B is noncompact in the q plane, it is often 
more convenient to describe the boundary in the complex z or y planes, where 
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z = - (1.12) 

and 

The variable y is commonly used in large-g series expansions. Some useful relations between these variables 
are 

9 = 1 + -, * = t^-> y = T^~ ( L14 ) 

v l+y l-z 

We define polar coordinates as 



z = (e ie (1.15) 



and 



y = pe lP (1.16) 



We shall use some standard notation from graph theory and combinatorics: v(G) and e(G) denote the 
numbers of vertices and edges in the graph G, and the symbol is given by 

qM=p\( q ) =T[{q-s) (1.17) 

The terms "edge" and "bond" will be used synonymously. For any partition of n, n = p + u, 

q ip) {q-p) [u) =q {u) {q-u) { rt = q^ +u ^ (1.18) 
For our results below, it will also be convenient to define the polynomial 

Dk{q) - = °" £(-«)-' = EHrf; 1 ) 9 ^ (L19) 

qKq ' j=0 s=0 ^ ' 

where 

a = q-l (1.20) 

and P(Ck, q) is the chromatic polynomial for the circuit graph, 

P(C k ,q) = a k + (~l) k a (1.21) 

Some useful properties of are listed in an appendix. 

The organization of the paper is as follows. In Section II we review some methods and results from our 
previous work on families of the form (-Kp)b + G r that will be necessary for our current discussion (where 
K p is the complete graph on p vertices p8[ and, following standard notation in the mathematical literature 
on graph theory, we use the symbol G + H for the "join" of these graphs; see below for definitions). Section 
III contains another method of constructing families of graphs with noncompact W boundaries B involving 
the removal of bonds not just from one vertex of the K p subgraph in K p + G r - In Section IV we present 
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and discuss a number of ways to construct families of graphs with noncompact W boundaries based on 
homeomorphic expansion of starting sets of graphs. In Section V we analyze two such methods involving 
homeomorphic expansion of the K p subgraph in a larger graph. Sections VI- VIII are detailed analyses of the 
respective chromatic polynomials, W functions, and their boundaries of regions of analyticity B for families 
of graphs obtained by specific homeomorphic expansions. In Section IX we include some remarks on general 
geometrical features of the families of graphs with noncompact W boundaries B. Finally, Section X contains 
our conclusions. 



II. FAMILIES OF THE FORM (K p ) b + G r 



Since we will construct and analyze a number of new families of graphs with noncompact boundaries B 
as homeomorphic expansions of families studied in Ref. ]lq ] it is first necessary to review briefly the method 
that we formulated and used in Ref. jl5| to generate such families. Consider a family of graphs G and its 
v(G) — > oo limit, {G} . If this family already has the property that the limiting function W({G},q) has a 
region boundary B that extends to complex infinity in the q plane, then we have no work to do to get such a 
boundary. So assume that this family is such that W({G}, q) has a region boundary B that does not extend 
to complex infinity in the q plane. Specifically, we start with an r-vertex graph G r and then form the graph 
K p + G ri where, as above, K p denotes the complete graph on p vertices ps| ], and we adopt the notation 
G + H that is commonly used in graph theory to indicate that every vertex of G is connected by a bond to 
every vertex of H; this is called the "join" of G and H in the mathematical literature. We then remove b 
bonds from one vertex in the K p subgraph; the resultant family is denoted 



(K p )b + G r 



(2.1) 



(In Ref. |Tf| we used the notation G x H for what is called G + H here, and the notation (K p x G n ) r b for 
the family denoted (K p )b + G r above, where the subscript rb signified "removed bonds".) Since each such 
vertex has degree A = p — 1 (i.e., has p — 1 bonds connecting to it), it follows that 



1 < b <p- 1 



(2.2) 



It was proved in Ref. ]l5j that in the limit r — > oo, the family (K p )b + G r has a noncompact W region 
boundary B(q), extending infinitely far from the origin of the q plane. A number of families of this type 
were constructed and the chromatic polynomials, resultant W functions, and boundaries B calculated [ fL5| . 
Specifically, by means of the relation 



P(K p + G r ,q) = q^P(G r ,q-p) 
the chromatic polynomial of (K p )b + G r was calculated to be |l5| 
P((K p ) b + G r ,q) = P{K P + G r , q) + bP{K p _ x + G r ,q) 



(2.3) 



= q^V 



Substituting the expression (1.10) yields 



P((K p ) b + G r ,q)=q^-V 



(? - (p - l))P(Gr, q-p) + bP(G r ,q - (p - 1)) 



iV 

(q - (p - l)){c (q - p) + g c 3 (q - p)aj (q - pf^ } 



(2.4) 
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+b{c (q - (p - 1)) + g Cj(q - (p - l))a,( g - (p - 1))*W} ! 2.:, ! 



For a given g, the boundary B as r — ► oo (with p fixed) is determined by the degeneracy in magnitude of 
the leading terms in eq. ( [2.5| ) at this value of q, viz., 

\a e (q-p)\ = \a e (q-(p-l))\ (2.6) 



This equation, and hence also B are independent of b for the interval (2.2), 1 < b < p — 1. In passing, we 
note that since the total number of vertices of (-Kp)b + G r is v((K p )i, + G r ) = p + r, another way to take 
n — > oo is to let p — > oo with 7' fixed, rather than letting 7' — > oo with p fixed. However, from the viewpoint 
of either statistical mechanics or graph theory, this is not as interesting a limit, since for any given graph 
G r and for any given (finite) value of q € Z + , as p becomes sufficiently large, one will not be able to color 
the graph (K p )b + G r and the chromatic polynomial will vanish. We shall therefore restrict ourselves to the 
other way of letting the number of vertices go to infinity, viz., r — > oo with p fixed. 
Writing 



at 



E > ( 2 - 7 ) 



recalling that the basic theorem that the coefficient of the highest-order term, q n , in the chromatic polynomial 
P(G n , q) of any n- vertex graph G n is unity implies that 

= 1 . (2-8) 



next dividing eq. (2.6) by |g ,Smax |, and finally reexpressing the degeneracy equation in the more convenient 



variable z given in eq. ( 1.12| ), we obtain 



|(l-pz) s — + J2 ««,.(! -J«)"*' 



|(l-(p-l)z) s — + £ a,, s (l-(p-l)z) s z s — ~ s \ (2.9) 



s=0 



This equation is clearly satisfied for z = 0, which shows that the boundary B for the r — > oo limit of the 
family of graphs {K p )i, + G r is noncompact in the q plane, passing through z = 1/q = 0. As noted above, 
it is thus more convenient to analyze B in the z plane. We next show that B passes vertically through the 



origin of the z plane for the r — > oo limit of all of the families of the type iK v )b + G r . Writing eq. (2.E) in 
polar coordinates, using ( |1.15| ), we have, for \z\ = £ — ► 0, 

(s max cos6 + 0(( 2 ) =0 (2.10) 

Since s max ^ 0, it follows that 

cos6»->±! as C^0 (2.11) 

which proves the above assertion, that B passes vertically through the point z = 0. 

The simplest families of graphs of the form (K p )t, + G r are those for which P(G r , q) has, in the notation 



of eq. ( 1.10 ), cq = 0, N a = 1, and ai(q) a linear function of q, 
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(2.12) 



For such families, in the q plane, B is a vertical line with arbitrary and 

g R = g c = V - + a ) for s moa; = 1 
where Re(q) = q R and Im(q) = q t . In the z plane, B is the circle defined by 



(2.13) 



(2.14) 



ith 



+ olq) 



(2.15) 



Equivalently, B is the circle \y — = ^ in the y plane, where y c = l/(p — (4 + ckq))- Thus, £> divides the 
z (equivalently q or y) plane into two regions: (i) R\, as defined above to include the region on the positive 
z axis contiguous with the origin, and (ii) its complement, denoted i?2- We find that 



W[ lim {{K p ) b + G r ],q) = q -p+ 1 + a for q G Ri 



(2.16) 



\W{lim[(K p ) b + G r ],q)\ = \q- P + a \ for q £ R 2 



(2.17) 



These families include the cases G r = (i) edgeless graphs K r ; (ii) tree graphs T r ; and (iii) chains of triangles 
with each pair of adjacent triangles intersecting on a mutual edge, i.e. {Ch)^, r in the notation of Ref. flBTl 
(The general chain graph of n p /c-sided polygons with adjacent fc-gons intersecting on a mutual edge is 
denoted (Ch)k, r ] the number of vertices, r, is given by r — (k — 2)n p + 2.) The chromatic polynomials for 
these families are 



P({K p ) b + K r , q) = (q - p) r + b(q - (p - 1)) 



(2.18) 



P((K p ) b + T r , q) = q^+V (q~(p+ 1))^ + b(q - pf 



(2.19) 



P((K p ) b + (Ch) k , r ,q) = q (p+1) \(q - (p + l))D fc (g ~ p) n * + bD k (q - (p - 



Thus 



(2.20) 



a = 0, -1, -2 for G r = K r , T r , {Ch) 3 , r 



(2.21) 



so that 



I i 3 

<? c = P - 2 ' P + 2 ' p + 2 fol Gr = ^ ' Tr ' ( C ' /l ) 3 ^ 



(2.22) 



The boundaries B for other families have the same behavior near the origin, z = 0, as shown by eq. (2.10), 
but have different features where they cross the positive real z axis. For example, consider G r = C r , which 
has, in the notation of eq. ( l.lCD , N a — 1 and a linear ai(q) = q — 1 but also a nonzero cq term, viz., 
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Co = (—l) r (q — 1). For lim r _» 00 [(/< p )fc + C r ], the portion of B near to z — forms part of the circle ( 2.14 ), 
but to the right of the points Zi nt . = 1/qint. and z* nt , where 

q in t.=P+2+— (2-23) 
the boundary bifurcates into two arcs, which cross the positive z axis at 

z c = — = — !— for G r = C r (2.24) 
q c P + l 

and at z = I /p. Since we shall compare this in detail with homeomorphic expansions of this family, we show 
this boundary in Fig. [|. For this and similar plots given later we also show chromatic zeros calculated for 
a reasonably large finite value of r (in this case, r = 27. The comparison of these chromatic zeros with the 
r —* oo locus B shows quantitatively how B forms via the merging of the zeros. Of course, since for finite r, 
the chromatic polynomial is of finite degree (equal to the number of vertices on the graph), its zeros (i.e., 
the chromatic zeros of the graph) are of bounded magnitude in q and hence do not track the asymptotic 
curve B all the way in to the origin of the z — 1/q plane. The same comment applies to the chromatic zeros 
that we will show for other families below. 




0.0 0.2 0.4 0.6 0.8 1.0 

Re(z) 



FIG. 1. Boundary B in the z = 1/q plane for the r — » oo limit of the family of graphs (ifa)b=l +C r = K2 + C r . Chromatic 
zeros for r = 27 are shown for comparison. 

Similar comments apply for other families such as (K p )b + L r and for the k > 4 members of the family 
(K p )i, + (Ch)k, r [ P"5| ; where L r denotes a cyclic square strip of width one square, i.e. cyclic ladder graph, 
and (Ch)k. r was defined above in our discussion of the special case k = 3. 
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Of the families of graphs that we have discussed, with noncompact boundaries B of regions of analyticity 
in W, most have chromatic polynomials that depend on three parameters, p, r, and 6, and the resultant B 
for the limit r — > oo depends on one parameter, p (and is independent of b). The family {K p )b + (Ch)k t r 
illustrates a case where the chromatic polynomial depends on four parameters - p, k, r, and b - and B thus 
depends on two parameters: p and k. 

Here we present results on another such family, 



G 



p,b,s,u 



(2.25) 



where sK u denotes the disjoint union of s copies of the K u graph (where the vertices of each of the s K u 
subgraphs are not connected to each other by bonds). That is, one adjoins K p to G r = sK u and then 
removes b bonds from one vertex of the K p subgraph. An illustration of the special case p = 2, b = 1, s = 3, 
and u = 2 is shown in Fig. @(d). 






FIG. 2. Illustrations ol lamilies of graphs considered in the text: (a) r = H EGk~2(K'2 + TV); (b) 
C fcir = HEG k _ 2 (K 2 + CV); (c) HEG k ^{K 3 + T r ); all for k = 3 and r = 4; (d) K 2 + sK u for s = 3 and u = 2. 



Clearly, the number of vertices in the disjoint union sK u is r = v{sK u ) = su, so that 

n = v((K p ) b + sK u ) = p + su (2.26) 

For s = 1 and 6 = 0, one has the simplification ifp + i^M = K p + U . Henceforth, we shall assume that b satisfies 
the restriction ( |2.2| ). For general s and u = 1, 

(^) 6 + s J ff 1 = (X p ) 6 +F s (2.27) 

so that our analysis for the latter family applies. For general u and s, from eq. ( |2.4| ) we have 

P((K p ) b + a*,,, g) - {(q-(p- l))[{q - p) (u) ] s + b[(q - (p - l)) (u) ] s } (2.28) 

This chromatic polynomial has the form 

P((K p ) b + sK u , q) = qV> \(q - p)^] V 6 , 5 (2.29) 
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where P bs is a polynomial in q of degree s. 

Now consider the limit s — > oo with p and u fixed. Equation ( [2.28 ) has the form of eq. ( l.ld| ), so that the 
degeneracy equation is 



|(g-p)M| = |( 9 -(p-l))H| 



Taking into account that q — is never a solution of eq. ( |2.30 ) and dividing by \q\ 
corresponding degeneracy equation of leading terms in the z plane: 

u— 1 u— 1 

i - (p+j»i - 1 n c 1 - (p- 1 

j=0 j'=0 



(2.30) 
we obtain the 

(2.31) 



By the same reasoning as before, it follows that the resultant locus B which is the solution to the degeneracy 
equation is noncompact in the q plane, passing through z = 1/q = 0. The general solution to the degeneracy 
equation ( [2.3C ) in the q plane is the vertical line given by 



Re(q) 

or equivalently, in the z plane, the circle 



p-l 



u 
2 



2 



Zc 

2 



p-l 



(2.32) 



(2.33) 



and similarly, in the y plane, the circle \y— | = tt where y c = l/(p— 2+~). The q (or z or y) plane is divided 
into two regions by the boundary B, namely, R\ with Re(q) > p — 1 + and R2, with Re(q) < p — 1 + ^. 
The dependence of q c on u is easy to understand: clearly, just as increasing p reduces the number of colors 
available to color the rest of the graph and hence increases the value of q c (serving as the demarcation on 
the real axis between the region of q values with W > 1 and nonzero ground state entropy Sq = ks In W, 
from the region where the behavior of W is different), so also increasing u for fixed p has the same effect. 
We find that 



W(lhn[(K p ) b + sK u \,q) = l(q-(p-l))M]V u for q 6 ifc 

s — *-oo 

where we recall the definition fll.l7| ). Further, 

\W(]im [(K p ) b + sK u },q)\ = \{q-p) {u) \ 1/u for q E R 2 

s — >oo 

It follows that 



W(lim[(K p ) b + sK u },q = q c )=0 if q c G 



if 



u is even 



(2.34) 



(2.35) 



(2.36) 



Note that the zeros in W extend into region to the right of q c , if p + u — 2 > q c , i.e, if u > 3. For 
example, for the case u — 3, for which q c = p + 5, the expression ( 2.34 ) has formal zeros atg=p — 1, q = p, 
and g = p + 1 . The first two of these are not relevant since they lie to the left of q c in the region R2 where 



the expression (2.34) does not apply; however, the third zero, at q = p + 1, is in region Ri and is a true 



zero of W . For this illustrative value u — 3, the expression (2.35) has formal zeros at q — p, q = p + 1, and 
q = p + 2. The last two of these are not relevant since they lie to the right of q c in region R\ where the 



expression (2.35) does not apply, while the first one, at q = p, lies in region Ri where this expression does 
apply and hence is a true zero of W in this region. 
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It is interesting to comment on the comparative values ol q c defined in the limit as r — > oo and the 
chromatic number x f° r specific members of a family of graphs parametrized by r. Before doing this, one 
should note that there are cases where B is trivial or does not cross the real axis, so that no q c is defined; 
an example is given by tree graphs T r , for which B = and \ = 2. For circuit graphs C r , x — 2 if r is 
even and x — 3 if r is odd. In the limit r — > oo, B is the circle \q — 1| = 1, so that g c = 2; hence this 
family illustrates the possibilities q c — X an d q c < X- For the square lattice with free boundary conditions or 
(periodic or toroidal) boundary conditions that preserve the bipartite nature of the lattice, x — 2 while in 
all cases q c — 3 |^3| , so that this illustrates the possibility q c > X- For the infinitely long strip of the square 
lattice with width L y = 3 vertices (w — L y — 1 squares), q c = 2, equal to x f° r this family; for the infinitely 
long strip of the triangular lattice with L y = 3, q c ~ 2.57, less than the chromatic number x = 3 for this 
family, while for the infinitely long triangular strip with width L y = 4, B does not cross the real q axis, so 
no q c is defined [|l3],[l8) . Clearly, there is a large variety of behavior as regards the question of whether q c for 
the infinite-vertex limit of a family of graphs is smaller than, equal to, or larger than the chromatic number 
X for specific members of the family. As an illustration of the families discussed here, K 2 + K r has x = 2 
and q c — 3/2 < \\ K 2 + T r has x = 3 and q c = 5/2 < x, and K 2 + (Ch) 3tr has x — 4 and q c = 7/2 < \- We 
shall discuss the situation for homeomorphic expansions of graphs in these families below. 



III. FAMILIES OF THE FORM {K p ) {b} + G r 

Another way to construct families with noncompact W boundaries B is to remove bonds not just from one 
vertex of the K p subgraph but also from one or more other vertices of this subgraph that are not adjacent 
to the first vertex, i.e., are not connected to this first vertex by any bonds. We symbolize this family by 

(K p ) {b} +G r (3.1) 

where the subscript {6} refers to the removal of multiple bonds from non-adjacent vertices of the K p subgraph. 
For this family, the maximum number of bonds that can be removed from the K p subgraph is the number 
of bonds in this subgraph, namely, 

e(K p ) = (t)= P -K^ (3-2) 



v 2y 

In general, the chromatic polynomial for such families is not quite as easily derived as that for the type of 



families where all bonds removed are from one vertex, where we obtained the simple formula (2.4) ]15| . 

A particularly simple family of graphs is obtained when one removes all of the (2) bonds in the K p 
subgraph, thereby yielding the subgraph K p and the resultant family 

K P + G r (3.3) 

Clearly, for p = 2, 

(K 2 ) b =i + G r = K 2 + G r (3.4) 

For the general family K p + G r , one can use the addition-contraction |3^] theorem to reexpress P(K p + G r , q) 
as the linear combination 
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P(K p + G r , q) = J2 S P )P ( R J + G n l) 



3=1 
P 



52sjpqWp{G r ,q-j) 

3=1 



(3.5) 



where the coefficients Sp^ are the Stirling numbers of the second kind, defined by the equation (cf. eq. 



0') 



(3.6) 



A closed-form expression is [Rll 



S U) 
p 



■EC- 1 )- 



fc=0 



fc 



fc" 



(3.7) 



Given that P(G,q) has the form (1.10), it follows, by the same argument as was used for the (K p )t, + G r 
families with removal of b bonds connecting to a given vertex, that B is noncompact in the q plane, passing 
through z = 1/q = 0. 

It is worthwhile to consider some special cases for G r . For G r = K r , we calculate 



P(K p +K r ,q)=J2s^q {j) (q-jY 

3=1 

Hence, in the limit r — > oo, B is determined by the degeneracy equation 

\q-l\ = \q-p\, i.e., |1 - z\ = |1 - pz\ 

Clearly, z = is a solution, and B is, in the q plane, the vertical line with 

p+l 



Re(q) = q c = 



i.e., in the z plane the circle 



2 



p+l 



(3.8) 



(3.9) 



(3.10) 



(3.11) 



Let us compare the values of q c for the respective r — > oo limits of K p + K r and (K p )b + K r (where it is 
understood in the latter case that b bonds are removed from a vertex in the K p subgraph). For p = 2, the 
value of q c in eq. ( 3.1CQ is equal to that for G r = K r in eq. (2.22), as follows from eq. (3.4). For p > 3, the 
value of q c in eq. (3.1C) is smaller than that for G r = K r in eq. ( 2.22Q . This is a consequence of the fact that 
the graph {K p )i, + K r can be obtained from K p + K r by addition of bonds and hence, by the theorem in 
section VII of Ref. |l7j , the coloring of the former graph with q colors is more constrained than the coloring 
of the latter graph; the same theorem implies that for q G Z + , 



W( lim [(K p ) b + K r ],q) < W( lim [K p + K r },q) 



(3.12) 



In region Ri to the right of the vertical line (3.10) in the q plane, 
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W( lim (K p + K r ),q) =q-l for qe Ri 

r—>oc 



(3.13) 



while in the complement, R2, 

\W(]im(K p + K r ),q)\=\q-p\ for q e R 2 

r— *oo 

A second case is G r — T r , i.e., the family K p + T r . We find 

P(K P + T r , q)=J2 4 3) 1 U) (l - - J - 1 

Hence, in the limit r — > 00, 2? is determined by the degeneracy equation 

|g-2| = \q- P -l\ , i.e., |1 - 2z\ = |1 - (p+ l)z\ 
In the (7 plane the locus of solutions to this equation is the vertical line with 

D ( \ P +3 

Re{q) = q c = — — 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



which again is a circle in the z plane, 



\z = — , z r = 

1 2 1 2 ' p+3 



(3.18) 



For p = 2, the value of <7 C in eq. ( 3.1 7|) is equal to that for G r — T r in eq. ( 2.22 ), as follows from eq. 
(3.4). For p > 3, q c is larger for linv^oo [(i4T p )t, + T r ] with b bonds removed from one vertex in K p than for 
linv^oo^p + T r ], for the reason given above. For region R% to the right of the vertical line (3.17) in the q 
plane, 



W( lim [K p +T r ],q)=q-2 for q e Ri 

r — >oo 



(3.19) 



while in the complement, i?2, 



\W{]im[K p + T r ],q)\ = \q-p-l\ for q & R 2 



(3.20) 



IV. CONSTRUCTION OF FAMILIES OF GRAPHS WITH NONCOMPACT B VIA 

HOMEOMORPHIC EXPANSION 

A major result of the present paper is the construction of families of graphs with noncompact W boundaries 
by means of homeomorphic expansion {HE) (also called inflation) of a beginning set of families of graphs. 
We recall the definition that two graphs G and H are homeomorphic to each other if H, say, can be obtained 
from G by successive insertions of degree- 2 vertices on bonds of G |ll]]. Each such insertion subdivides an 
existing edge of G into two, connected by the inserted degree-2 vertex. This process is called homeomorphic 
expansion and its inverse is called homeomorphic reduction, i.e. the successive removal of vertices of degree 
2 from a graph H. Clearly, homeomorphic expansion of a graph always yields another graph. The inverse 
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is not necessarily true; i.e., homeomorphic reduction of a graph can produce a multigraph or pseudograph 
instead of a (proper) graph 0. Here, a multigraph is a finite set of vertices and bonds that, like a graph, 
has no bonds that loop around from a given vertex back to itself but, in contrast to a (proper) graph, may 
have more than one bond connecting two vertices. A pseudograph is a finite set of vertices and bonds that 
may have multiple bonds connecting two vertices and may also have looping bonds. For example, consider 
homeomorphic reduction of a circuit graph C' r : removing one of the vertices (all of which have degree 2), one 
goes from C r to C r _i, and so forth, until one gets to C3. During this sequence of homeomorphic reductions, 
one remains within the category of graphs. However, the next homeomorphic reduction takes C3 to C2, 
which is a multigraph, not a proper graph. Removing one of the two vertices in C2 produces a pseudograph 
consisting of a single vertex and a bond that goes out and loops back to this vertex. Thus homcomorphism 
is an equivalence relation on pseudographs. This complication will not be relevant here because we shall 
only use homeomorphic expansions, not reductions, of graphs, and the homeomorphic expansion of a proper 
graph always yields another proper graph. For our subsequent discussion, we shall denote the homeomorphic 
expansion involving the insertion of fc— 2 additional vertices (where fc > 3) on a specific bond & of a graph G as 
HEGk-2\b{G). Most of our studies will be of graphs in which the homeomorphic expansion of G is performed 
systematically on each bond of G; in these cases, we shall denote the resultant graph as HEGk-i{G). 

The homeomorphic expansion of a family of graphs with a compact (empty or nontrivial) locus B can 
produce a family of graphs with either a compact or noncompact B; this depends on the nature of the original 
family and of the homeomorphic expansion. For example, if one starts with the tree graph T r and inserts 
fc — 2 degree-2 vertice on each of the bonds, one obtains another tree graph, T r > , where r' = r + (r — l)(fe — 2). 
In the limit r — ► 00, both T r and its homeomorphic expansion, T r /, have a trivial B = 0. If one starts with 
the circuit graph C r and inserts fc — 2 degree-2 vertices on each bond of this graph, one obtains another 
circuit graph, C r i, where r' = r + r(fc — 2). In the limit r — > 00, both C r and its homeomorphic expansion 
C r ' have a compact B given by the unit circle centered at # = 1, \q — 1| = 1. We next proceed to discuss the 
cases of main interest here, where the homeomorphic expansion (i) leads from a family with a compact B to 
one with a noncompact B or (ii) takes a family that already has a noncompact B to another that again has 
a noncompact B. We comment on the effect of a homeomorphic expansion of a graph G on its girth 7(G), 
defined as the length of (= number of vertices on) a minimal-length closed circuit in this graph. Clearly, a 
homeomorphic expansion applied to one or more bonds of a graph increases the girth of the graph if and 
only if these bonds lie on the minimal-length circuits of the original graph. We shall consider four main 
types of homeomorphic expansions, as well as combinations thereof: 



1. HEKO 



Start with the family K p + G r (where no bonds have been removed from any vertex of the K p subgraph) . 
If in the limit r — > 00, B for the G r family itself is compact (bounded) in the q plane, then the same holds 
for the family K p + G r . Now insert degree-2 vertices on one or more bonds of the K p . We shall denote this 
generically as an HEKO homeomorphic expansion, meaning that the homeomorphic expansion acts on the 
K p subgraph and that originally there were 6 = bonds removed from this K p . Specifically, we denote the 
graph obtained by successive insertion of fc — 2 degree-2 vertices on a single bond bj of the K p as 

HEK k _ 2;b .(K p + G r ) (4.1) 
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where k > 3, and so forth for similar homeomorphic insertions on other bonds of the K p . We will show 
below that this homeomorphic expansion leads from the family K p + G r with compact (trivial or 
nontrivial) B to the family fl4.ip with a locus B that is noncompact in the q plane. The labelling convention 
is chosen so that as one moves along the expanded set of bonds linking a vertex of the original K p to what 
was originally an adjacent vertex, one traverses a total of k vertices, including the original pair, i.e., k — 2 
inserted vertices. 



Next, we consider cases (nos. 2-5 below) where we begin with a family of graphs that already has a locus B 
that is unbounded in the q plane. For these the homeomorphic expansion produces another family of 
graphs again with an unbounded B: 

2. HEK 

Start with the family {K p ) b + G r where b bonds have been removed from one vertex of the K p subgraph 
and insert degree-2 vertices on one or more bonds of the K p . We shall denote this generically as an HEK 
homeomorphic expansion, meaning that the homeomorphic expansion acts on the K p subgraph. 
Specifically, we denote the graph obtained by successive insertion of k — 2 degree-2 vertices on a single 
bond of the K p as 

HEK k ^ 2 , b] [(K p ) b + G r ] (4.2) 
where again k > 3, and so forth for similar insertions on other bonds of the K p subgraph. 

3. HEG 

Start with either (K p ) b + G r or (K p )^ + G r and add vertices to bonds in the G r subgraph. We label this 
type of homeomorphic expansion generically as HEG. Analogously to the previous category, we denote the 
respective graphs obtained by successive insertion of k — 2 degree-2 vertices on a single bond of the G r 
subgraph as 

HEG k - 2 ; bj [{K p ) b + G r ] , HEG k - 2 ;b 3 [{K p ) {b} + G r ] (4.3) 
and the graphs obtained by successive insertions of k — 2 degree-2 vertices on all of the bonds of G r as 

HEG k - 2 [{K p ) b + G r ] , HEG k - 2 [(K p ) {b} + G r ] (4.4) 



4. HEG 

Start with either (K p ) b + G r or (K p ){ b y + G r and add vertices to the bonds connecting vertices in the K p 
subgraph to vertices in the G r subgraph. We label this type of homeomorphic expansion generically as 
HEC, where the "C" refers to the fact that the homeomorphic expansion is performed on the 
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above-mentioned connecting bonds. We denote the respective graphs obtained by successive insertion of 
k — 2 degree-2 vertices on a single bond bij connecting a vertex Vi £ K p to a vertex Vj € G r as 

HEC k - 2;bij {(K p ) b + G r ] , HEC k >:/, + Gr] (4.5) 

and the graphs obtained by successive insertions of k — 2 degree-2 vertices on all of the bonds connecting 
vertices of K p to vertices of G r as 

ff£C^ 2 [(X p ) b + G r ] , ff£C fc _ 2 [(tf p ) {6} + G r ] (4.6) 



5. Combinations 



Clearly, one can combine several types of homeomorphic expansion. For example, starting with 

(K p )b + G r , one can add vertices both to bonds in the subgraph K p , to bonds that connect K p to G r , and 

to bonds in the G r subgraph. 

In the present paper we shall concentrate on homeomorphic expansions of types (l)-(3). Our results for 
homeomorphic expansions of type (4) involve somewhat more complicated boundaries B than those discussed 
here and will be dealt with in a separate paper. Composite homeomorphic expansions of type (5) can be 
studied by similar means. We now proceed to consider the various homeomorphic classes in more detail. 

Given that the number of vertices in the original graph is a linear function of the two (positive integer) 
parameters p and r, the number n of vertices of the homeomorphic expansion is a linear function of p, r, 
and k. There are therefore three basic ways of producing the limit n — > oo, namely (L denotes limit) 

L p : p — > oo with r and k fixed (4-7) 



L, 



oo with p and k fixed 



(4.8) 



L k : k 



with 



P 



and 



fixed 



(4.9) 



We have explained above (after eq. (2.6)) why the limit L p is not very interesting. From the viewpoint of the 
present work on boundaries B that are noncompact in the q plane, the limit is also not of primary interest, 
since it generically yields a compact boundary B, as we shall illustrate below. Hence we shall concentrate 
on the limit L r here. 



V. HEKO AND HEK HOMEOMORPHIC CLASSES 

Let us consider first the homeomorphic expansions of type (1), namely, HEKk-2.bj (Kp + G r ). By use of 



the deletion-contraction theorem and eq. (2.3), we find that, for arbitrary G r , 



P[HEK k ^ bj (K p + G r ),q) = D k P(K p + G r ,q) + [D k + (-l) fc_1 ]P(A" p _i + G r , q) 

= D k q^P(G r ,q-p) + [D k + (-l) k - 1 ] q ^P(G r ,q-(p-l)) (5.1) 
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It follows by the same argument as that given with eq. that as r — > oo, this family HEKk-2>i,b4 (Kp + 

G r ) has a W boundary B that is noncompact in the q plane. This is true independent of whether the family 



G r has this property. Indeed, as is evident from a comparison of eq. (5.1) with eq. (|2.4[), for a given G r , 



B( lim HEK k - 2 ,bAK p + G r j] = b( Urn [(K p ) b + G r }) (5.2) 

That is, if we start with the family K p + G r , remove b bonds from a given vertex in the K p subgraph, and 
take r — ► oo, the resultant boundary B is the same as if we had instead homeomorphically added some 
number k — 2 > 1 degree-2 vertices to a bond in the K p subgraph. One can, of course, continue this process 



with homeomorphic expansions of other bonds of the K p . Note that if one uses the limit L k in eq. (4.9) to 



get n — > oo for this class of homeomorphic expansions, then, for arbitrary (finite) G ri the resultant boundary 



is the compact locus comprised by the unit circle \a\ = 1 (where a was defined in eq. (1.20), i.c 



2?L fc = q such that \q — 1| = 1 (5.3) 



We next consider the homeomorphic expansions of type (2), namely, HEKk-2;bi i(K p ) b + G r ]. Here, the 
boundary B for the r — > oo of the beginning family (K p )b + G r is already noncompact in the q plane, and 
the families generated by the homeomorphic expansion maintain this property. As an illustration, consider 
the family (K^)i + G r and homeomorphically expand one of the two remaining bonds, denoted bj, in the 
(^3)1 = T3 subgraph. We find for the resultant chromatic polynomial the result 

P(HEK k ^ 2 , b] [(#3)1 + G r ],q) = D k q&P(G r ,q - 3) 

+ [2D k + {-if-^PiGr, q - 2) (5.4) 

Again, using the same reasoning that we employed before with eq. (|]J), we deduce that as r — > 00, this 
family HEKk-2;bj[(^-3)i + Gr] nas a noncompact locus B in the q plane. As was the case for the HEKQ 
class, for the limit L k , the resultant boundary is the unit circle \q — 1| = 1 given in eq. ( |5.3| ). 



VI. THE FAMILY T fc , r = HEG k - 2 (K 2 + T r ) 



We next proceed to analyze in detail several homeomorphic classes of graph families of the form HEG. 
One interesting family of graphs with noncompact B{q) is obtained by homeomorphic expansion starting 
with the family (K p )b + T r for r > 2 and adding vertices on each of the bonds in the T r subgraph. We let k 
be the number of vertices on T r between each pair of vertices that connect with the K p , including this pair 
(which were originally adjacent on T r before the homeomorphic expansion). We denote this family as 

T p , b , k ,r = HEG k - 2 [(K p )b + T r ] (6.1) 

where as above, r > 2 and 1 < b < p — 1. The number of vertices is given by 

v{T pA k, r ) = (r - l)(fe - 1) + p + 1 (6.2) 

For the lowest value k = 2, 

T p ,b,k=2,r = (K p )b + T r (6.3) 
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which we studied previously fl5| . We shall thus concentrate on the cases k > 3 here. It suffices for our 
present purposes to consider the simplest nontrivial case p — 2 and hence b = 1, for which 



(K 2 ) b=1 +T r = K 2 + T r 
An illustration of a graph of this type is given in Fig. |2|(a). For brevity, we define 

Tk,r — T p= 2fi=l,k,r = H EGk-2 (K 2 + T r ) 

We observe that the chromatic number is 

x(Tk,r) = 2 for k odd 
= 3 for k even 

The girth, i.e., the length of a minimal- length circuit on this graph, is 

l(Tk,r) =k+l 

By the deletion-contraction theorem, we find the recursion relation 

P(T fe , r , q) = [D 3 D k + (— l) fe_1 J P(T fc , r _i, ? ) + q(q - l)D k (D k+1 ) r 
Solving this, we obtain the chromatic polynomial 

P(T ktr ,q) = A k [D a D k + (-l)*" 1 ] ^ 2 + q{q - ^(Dfe+xf- 1 

where 

A k = (q-3)P(C k+1 ,q) + (q-l)P(C k - 1 ,q) 



(6.4) 



(6.5) 



(6.6) 



(6.7) 



(6.8) 



(6.9) 



(6.10) 



P(T kyr ,q) has the general factors q(q — 1) for k odd and q(q — l)(q — 2) for k even. As is evident from eq. 

(6.11) 



flOD, P(T ktr ,q) has the form of eq. Ql.lO| ) with 

«i = -Dfc+i 
and 

a 2 - ^3^fc + (-l) fe_1 



(6.12) 



(recall that D% = q — 2). In the limit r — > oo (with fe fixed [p2| ) the locus B is determined by the degeneracy 
of magnitudes 

|oi| = |a 2 | (6.13) 

This equation can be simply expressed in terms of the variable a = q — 1 which was defined above in eq. 
( 1.20 ) (and should not be confused with the variables ai(q) and 02(9)). To do this, we multiply both sides 
of eq. (ETl3|) by \q(q - 1)| = |(a + l)o| (the spurious solutions at 5 = and q = 1 thereby introduced are 
discarded) and use eq. ( 1.2l| ) to get 



\a(a k - a*- 1 + 2(-l) fe+1 )| = \a(a k + (-l) k+1 ) 



(6.14) 
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i.e., after dividing by \a 



fc+i 



ii- y +2(-i) fc +yi = |i + (-ir+vi 



ft+l„,ft| 



(6.15) 



Since y = is a solution of this equation, B is noncompact in the q plane, passing through y = z = 0. In 
polar coordinates, with y = pe l13 , eq. ( 3.15] ) yields 



3/r 



2fc-l 



+ p - 2cos/3 + 2(-l) fe -V i cos(fc/3) - 2 / 9cos((fc - l)/3) 



= 



(6.16) 



As p — > 0, it follows that cos/3 = 0, i.e. /3 = ±7r/2, so that i3 approaches y = z = vertically. 

To calculate the point at which the boundary B crosses the real q, or equivalently, y or z axes, we set 
[3 = in eq. (6.16); for p ^ this gives 



•V 



2 ft- 1 



2(-l) 



2p) = 



(6.17) 



(Note that Re(q) > Re(z) > and Re(q) > 1 Re(y) > so that setting /3 = 0, i.e. making ?/ real 
and positive, implies that q is real and q > 1; we comment later on the interval < <? < 1.) Since p = \y\ is 



real and positive, we are only interested in roots of eq. ( 6.17 ) that have this property. For both even and 
odd k, there is only one such root, which is thus y c (equivalently, z c ), the minimum nonzero value of y (z) 
respectively, at which the boundary B crosses the positive real axis in the y (z) plane, corresponding to the 
maximal finite real point q c at which the pre-image of this boundary B crosses the positive real axis in the q 
plane. For k odd, this single root is p = y = 1 (the expression on the left-hand side contains a factor (p— 1)) 
so that z c = 1/2, and 



2 = x f° r k odd 



(6.18) 



For k even, the single such root of eq. (3.17) has a value that increases monotonically from p — y c = 2/3 
(equiv. z c — 2/5) for k = 2 toward the limit p = y c = 1 (equiv. z c = 1/2) as k — > oo through even integers. 
Correspondingly, for k even, q c decreases monotonically from 5/2 for k = 2 toward 2 in the same limit; some 
values are q c = 2.2564, 2.1736, 2.1315 for k = 4, 6, 8). Hence, 



2 < q c < - < X = 3 for 



even 



(6.19) 



To show that B does not cross the positive real q axis in the interval < q < 1, we take [3 — n, i.e., y real 
and negative, in eq. ( |6.16 ), which, for p ^ 0, yields the equation 



3p 



2fc-l 



p + 2 - 2p k ~ 1 - Ap k = 



(6.20) 



This range — oo < y < includes both the interval < q < 1 and the interval — oo < q < 0. We find that in 



this range of y, eq. (6.20) has no real positive roots for p except for p = 1, i.e., y = —1 or equivalently, q = 0, 
but this is just a spurious solution introduced when we multiplied both sides of eq. ( |6.13j ) by \q(q — 1)| to 



express it as eq. (6.14). Hence, B crosses the positive z or equivalently q axis only once for this family of 
graphs. 

We proceed to discuss the boundary B for specific values of k further. For k = 2 (l5) B is the vertical line 
with Re(q) — 5/2 in the q plane, or equivalently, the circles \z — z c /2\ — z c /2, \y — y c /2\ = y c /2 in the z and 
y planes, where z c = 2/5, y c = 2/3 as given above. 
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FIG. 4. As in Fig. ^ for k = (a) 5 (b) 6. Chromatic zeros for r = (a) 10 (b) 8 are shown for comparison. 

The boundary B is shown in the z plane for 3 < k < 6 in Figs. [| and || [j33f . In region Ri containing the 
portion of the positive z axis analytically connected to the point z = 0, 



W{[]im T ktT ],q) = (oi) 



i/(fc-i) _ 



(^ fe+1 ) 1 /( fe - 1 ) for q&R x 



(6.21) 
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In region i?2 occupying the rest of the z plane for 2 < k < 4, and the rest of the z plane except for the 
interiors of the other closed loops for k > 5, 

\W([]imT k , r ),q)\ = \a 2 \ 1 /^ = \D 3 D k + (-l) k - 1 \ 1 ^ k - 1 ^ for q G R 2 (6.22) 

r— »oo 

For odd A;, as q decreases toward 2 one has, using eqs. ( ET2ll ) and ( pTi| ), 



lim W([ lim T fcr ],g) = 1 (6.23) 

<j\2 v — >oo 



and, as q increases toward 2, one has, using eq. ( |6.22|) , 



lim||W([lim Tfc, r ],g)| = l (6.24) 

q/*2 r— >oo 



The equality of the right-hand sides of eqs. (6.23) and (6.24) follows from the fact that q — 2 is on the 



boundary B for odd k. In each of the additional regions enclosed by closed curves, 

\W([]im T k , r ],q)\ = M 1 ^-!) = l^+il 1 /^-!) (6.25) 

r — >oo 

For comparison, we show chromatic zeros for reasonably large finite values r in each of these figures. 

We find several general features of the boundary B for the r — > oo limit of this family of graphs. The first 
feature is that B has support only for non-negative Re(z) or equivalently, non-negative Re(q): 

zeB^Re(z)>0 i.e., Re(q) > (6.26) 

Indeed, the only place on the boundary B where Re(z) vanishes is at the origin z = itself: 

{zeB and Re(z) = 0) => z = (6.27) 

Thus, in the vicinity of the point z = Q, the curve comprising B bends to the right as one increases Im(z) 
above zero or decreases it below zero; we have shown above that this curve crosses the point z — vertically. 
The second feature is that B has no multiple points ]34] ] . The third feature concerns the number of different 
connected components N comp . comprising B. In the present case, this is simply related to the number of 
regions N reg , by the equation N reg , = N comp ,+l. One might hope that there would be a general mathematical 
theorem that would state the number of different connected components N comp . of the solution set of a given 
polynomial equation without requiring an explicit solution. However this question, which is related to the 
sixteenth Hilbert Problem, still remains unanswered p6] , p7| . One can make use of an upper bound on -/V comp 
contained in the Harnack theorem, which is N comp < g+ 1, where g denotes the genus of the algebraic curve 
p5[ . For plane algebraic curves of the type relevant here, 

E«X = (6.28) 

m,n 

with maximal degree d, where y R = Re(y), y r = Im(y), the genus is g — (d — l)(d — 2)/2. Thus, for a 
given case, one writes the degeneracy equation ( |6.15 ) or its equivalent in the z plane, out into cartesian 



components. The case k = 2 has already been discussed and leads to the equation (y R — l/3) 2 + y^ = (1/3) 2 , 
or equivalently, (z R — 1/5) 2 + z 2 — (1/5) 2 where z R = Re(z) and z 2 = Im(z). These equations have 
homogeneous degree |55) d — 2, hence genus g = 0, so Harnack's theorem yields the bound N comp , < 1 
which, together with the fact that N comp , > 1 implies that N comp . = 1, in agreement with the elementary 
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solution of the explicit equation above. However, for higher values of k, the Harnack upper bound is not 
very restrictive. For k > 3, eq. ( |6 .15 ), when written out in cartesian components, is 



(y R + iy + yi\F(y R ,y 1 )=Q (6.29) 
where F(y R ,y I ) = 3(y R k ~ 2 + y 2k ~ 2 )+ lower order terms. Because the first factor is positive definite, eq. 



(S.15) thus reduces to F(y R ,y 1 ) = 0, of homogeneous degree d = 2k — 2 and hence genus g — (2k — 3)(fc — 2). 



The Harnack theorem then yields the upper bound N comp . < 2k 2 — 7k + 7. As an example, for k = 3, 4, 5, 6, 
this bound has the respective values 4,11,22,37, while the actual values are N comp , = 1,1,3,3, and so forth 
for higher k. 

For both of the cases k = 3,4, as was true for k = 2, B divides the z plane into two regions: (i) the 
region R\ including the vicinity of the positive real axis contiguous with the origin, z = 0; and (ii) the region 
i?2 occupying the rest of the z plane outside of R±. In the z plane, the chromatic zeros tend to cluster 
on the regions of the curve B in the "northeast" and "southeast" directions. Indeed, as k increases from 
3 to 4, B(z) develops protuberances in these "northeast" and "southeast" directions. Since for a given r 
and k the chromatic zeros are bounded, they avoid the origin of the z plane, as can be seen in the figures. 
Between k — 4 and k = 5 there occurs a qualitative change in B, namely, that, whereas it consisted of a 
single component for 2 < k < 4, it consists of three disconnected components for 5 < k < 7, comprised of a 
self-conjugate closed curve together with a complex-conjugate pair of closed curves. Hence, while B divides 
the z plane into two regions for 2 < k < 4, it divides this plane into four regions for 5 < k < 7. In a figurative 
language, one can think of the northeast and southeast bulges that are evident in the k = 4 case as breaking 
off to form the two separate closed curves that first appear for k = 5. We also observe that the additional 
closed curves and associated disconnected regions appear to line up approximately along the vertical line 
with Re(z) = 1/2. It is interesting to note that this line maps to unit circles in the y and q planes: 

z : Re(z) — —. Im(z) arbitrary ■<=>■ \y\ = 1 <^=^ \q — 1| = 1 (6.30) 

Thus, the disconnected phases that appear in the z plane are clustered around the unit circle \y\ = 1 in the 
y plane. In the q plane, portions of the boundary B are also clustered around the circle \q— 1| = 1 while one 
portion extends infinitely far from the origin. 



VII. THE FAMILY C k , r = HEG k - 2 (K 2 + CV) 

We obtain a further infinite family of graphs with noncompact B(q) by the same steps as in the previous 
section, but replacing the subgraph T r by the circuit subgraph C r . That is, we homeomorphically expand 
the family {K p \ + C r for r > 2 by adding vertices on each of the bonds in the C r subgraph. As before, 
we let k be the number of vertices on C r between each pair of vertices that connect with the K p subgraph, 
including this pair. We denote this family as 

C p , b , k ,r = HEGk-%[{K p ) b + C r ] (7.1) 

where as above, r > 2 and 1 < b < p — 1. For a given set of parameters p, k, and r, the number of vertices 
in this family is 

v(C pAkir )=r{k-l)+p (7.2) 
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For the lowest value, k = 2, 



C Pt b,k=2,r = (Kp)b + C r 



(7.3) 



which we studied previously accordingly, we concentrate here on the cases k > 3. It suffices for our 
present purposes to consider the simplest nontrivial case p — 2 and hence 6=1, for which 



{K 2 ) b =i + C r = K 2 + C r 

As before, a short notation is useful: 

Ck,r = Cp=2,5=l,fc,r = HEG k -2(K 2 + C r 



(7.4) 



(7.5) 



An illustration of a graph of this type is given in Fig. 0(b). The chromatic number and girth are the same 
as for Tk, r , given by the right-hand sides of eqs. ( p$) and (pTij ). By the deletion-contraction theorem, we 
find the recursion relation 



P(C k , r ,q) = D k P(T k , r ,q) + {-l) k - l P{C k , r . x ,q) 
For the lowest value of r, namely, r — 2, we have 

P(C k<2 ,q) = \D k + (-l) fc - 1 lp(T fci a,g) + (-1)'^)'^,?) 



(7.6) 



(7.7) 



Solving the recursion relation ( |7.6[ ) with ( |7.7[ ), we calculate the chromatic polynomial. As before, we express 
this in the form of eq. ( 1.10| ). (In order to render the polynomial property manifest, one divides through by 
the factors in the denominator, thereby generating a series, using the identity (x m - l)/(x - 1) = £™ o xK) 
We get 



P{C kir ,q) = D k 



A fc a 2 (a 2 + (-l) fc - 1 ) 



gfe-lWffli + Hf- 1 ) 



+ (-i) 



(fc-l)r 



P(Ck,2,q) 



(7. 



where Afe was defined in eq. ( 6.1 C| ) . P{C k}T , q) has the general factors q(q — 1) for k odd and q(q — l)(q — 2) 
for k even. P(C k>r , q) has the form of eq. ( 1.1 0| ) with a\ and a 2 as given in eqs. (6.11) and ( 6.1 2| ) , together 
with 



a 3 — 1 



(7.9) 



For large \q\ (small \z\), the terms a\ and a 2 have larger magnitudes than 03, so that in the limit r — * 00 (with 
fc fixed |^2| ) the boundary i3 is determined by the degeneracy equation |oi| = |a 2 |, and the same analysis 
goes through as before for the family HEG k -2(K2 + T r ), with the same conclusions that B is noncompact 
in the q plane, passing through the point z = y = vertically. For the lowest case k = 2, i.e., K2 + C r , 
B consists of three regions, as shown in Fig. |l|): R\ including the real axis in the interval < z < z c with 
z c = 1/3 as given by eq. ( [2.24| ); R 3 centered around z = 2/5 and occupying the interval z c < z < 1/2 on 
the real axis; and i? 2 occupying the complement of the z plane extending outward to \z\ — > 00. This region 
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diagram may be contrasted with the simpler one for K2 + T r , which was just a circle, separating the z plane 
into two regions. 

In Figs. |5| and ^| we show the boundary B in the z plane for 3 < k < 6 j33|. Again, for comparison, wc 
show chromatic zeros for reasonably large finite values r in each of these figures. The property that B passes 
through z — vertically is evident in these figures, as are the properties (6.26) and ( |6.27 ). For each value 



of k, the region R\ is the one occupying the non-negative interval of the z axis adjacent to the origin. The 
region R 2 occupies the rest of the z plane for 2 < k < 4, and the rest of the z plane except for the interiors 
of the other closed loops for k > 5. For odd k, there are two regions on the real z axis: R\ for < z < 1/2 
and i?2 for z < and z > 1/2, while for even k > 4 there are these phases and, in addition, a third "pocket" 
phase i?3 whose right-hand boundary with R2 occurs at z = 1/2 and whose left-hand boundary with R\ 
occurs at the point z — z c slightly less than 1/2 (see further below). 
As was the case for \im r ^ 00 [HEGk-2(K2 + T r )), we find 

WQlim C fc>r ],g) = (a 1 ) 1 /(fe- 1 ) =( J D fe+1 )V(fe-i) for q e Rl ( 7 . 10 ) 

r^oo 

and 

\W([lim C k , r ],q)\ = \a2\ 1/(k - 1) = \D 3 D k + (-l) k ~ 1 \ 1 /( k -V for q e R 2 (7.11) 



Hence for k = 3 the equalities ( 3.23] ) and (6.24) hold, with T k , r replaced by C k>r . For k > 3, additional 



regions appear as "outgrowths" on the right-hand boundary of R\ in the z plane, bounded on the left by 
Ri and on the right by i?2- We label these regions generically as Ri-2-j, where j > 3 indexes the region. 
For example, there are two such regions for k = 3, which are complex-conjugates of each other, and there 
are three such regions for k = 4, consisting of a complex-conjugate pair and the region i? 3 that includes 
an interval of the real z axis between z — 0.42495 and z = 1/2. For arbitrary k > 3, in these outgrowth 



regions lying between Ri and R2, the leading term in the limit r — » 00 is the last one in eq. (7.5), viz., 
(_l)(fc-i)-p(C fe!2 , (Z ), so that 

\W{[]im C k;r ],q)\ =1 for q G R\-2,j (7.12) 

r — >oo 

For k > 5, additional sets of regions occur that are not connected with Ri ; these consist of complex-conjugate 
pairs, lying near to the vertical line in the z plane with Re(z) — 1/2. For the cases that we have studied, 
each of these disconnected sets of regions actually consists of two: for the ones with Im(z) > 0, one part 
occupying a "northeast" position and the other a "southwest" position, as is evident in the figures. We label 
these regions disconnected from Ri as RdiscNE-j, Rdisc.swj, and their complex-conjugates as R* disc j^E-ji 
SW-j ■ We- find that 

\W([]im C fc , r ], (Z )| = |a 1 | 1 /( fe - 1 ) = | J D fe+1 | 1 /( fc - 1 ) for q e R disc , sw .j, R* Aisc , S w,j (7-13) 



\W{{ lim C k ,r], q)\ = 1 for q G R di8C , N E;j, R* disc NE; j (7-14) 

r — >oc ' ,J 

As was noted before in the case of the r — > 00 limit of the family T kt r, this clustering of the disconnected 
phases along the vertical line Re{z) — 1/2 is equivalent to their clustering around the unit circle in the y 
plane. 

For odd k, B crosses the real z axis away from z = (cquivalently the real q axis) once, at z — 1/2 (i.e., 
q = 2). For even fc, B crosses the real z axis at two points away from the origin z — 0: (i) at z = 1/2, and 
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(ii) at a value of z that increases monotonically from z = 1/3 for k = 2, approaching z = l/2asfc^oo 
through even integers. Hence, taking into account that 



we have 



while 



x(Ck,r) = 2 for k odd 
= 3 for k even 



q c = 2 = x f° r °dd k 



(7.15) 



(7.16) 



2 < q c _ 3 = x f° r even fc 



(7.17) 



Equations ( 7.16 ) and ( 7.17 ) may be compared with the analogous eqs. ( 3.18 ) and ( 6.1 9| ) for the family 
Tk : r = HEGk-2{K2 + T r ). To derive the above results, we observe first that for the case of odd fc, this 
follows from the same analysis as applied to the r — > oo limit of the family HEGk-2(K2 + T r ) in the previous 
section, since the crossing is determined by the degeneracy condition |ai| = |«2 1 - For even k 7 the crossing 
point at the smaller value of z is determined by the degeneracy condition of leading terms 



S{Ri,R 3 ) 



Kl = |«3 1 = i 



(7.18) 



To put this into a form that can be used for the analysis of our region diagrams in the z plane, Figs. ^| and 



|6j, we multiply both sides of eq. (|7.18| ) by \q(q — 1)| = \(a + l)a\ (the spurious solutions at q = and q = 1 
thereby introduced are understood to be discarded) . The resulting equation is most simply expressed in the 
y plane, as |1 - y k \ = \y k ^ 1 {y + 1)|, i.e., 



B(R U R 3 ): p 2k ~ 2 + 2p k (p k - 1 cosf3 + cos(k[3)j -1 = 



(7.19) 



For P = 0, this equation has a single acceptable (real positive) root for p = y c that increases monotonically 
from 1/2 (i.e., z c = 1/3) for k = 2 through 0.738984 (z = 0.424951) for k = 4, toward p = y c = 1 (z c = 1/2) 
as k goes to infinity through even values; cquivalcntly, q c decreases from 3 at k — 2 toward 2 as k — > oo 
through even values. Some explicit values are q c — 2.35321, 2.21486, 2.15442 for k = 4,6,8. The crossing 
point at the larger value of z (for even k), viz., z = 1/2, is determined by the degeneracy condition of 
leading terms \ai\ — \a 3 \. Multiplying by \q(q — 1)| and reexpressing this in terms of y as before, we obtain 



|l-2/ + 2(-l) fc - 1 y fc | = |(l+ 2/ ) 2/ fe - 1 |, 



i.e. 



B(R 2 ,R 3 ) 



3p 2k _ p 2k-2 



1 - 2p{p 2k - 2 + l)cos/3 



-Ap k 



cos(k(3) — pcos((fc — 1)(3) 







(7.20) 



For j3 = 0, this equation always has a root at p = y = 1, (i.e., z = 1/2, q = 2). (It also has a real positive 
root at a value of p < 1 (i.e., z < 1/2), but 02 is not a leading term in this region, so that this root is 
irrelevant.) 

The number of multiple points on B is 



N„ 



2(k — 1) , for even k 



(7.21) 
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comprised of k — 1 complex-conjugate pairs, and 



N m . p . = 2k - 3 , for odd k (7.22) 

consisting of a real one at z = 1/2, i.e., q = 2, together with fc — 2 complex-conjugate pairs. 

It is interesting to compare the boundaries B and associated region diagrams for the r — > oo limits of the 
two families of graphs HEGk-2(K 2 + T r ) and HEGk-2(K 2 + C r ) for a given value of k. We have done this 
above for k = 2. For fc > 3, one sees that the diagrams for HEGk-2(K2 + C r ) look somewhat similar to those 
for HEGk-2(K2+T r ) with the differences that (i) there are additional "outgrowth" regions contiguous to the 
right-hand side of R\ ; (ii) the additional disconnected regions that appear for k > 5 are themselves composed 
of additional regions; (iii) B contains multiple points; and (iv) for even k there are three, rather than just 
two regions along the real z axis. Owing to the presence of the multiple points on B (i.e. singular points, in 
the sense of algebraic geometry), the analysis of bounds on the number of disconnected components, N comp , 
is more difficult than for algebraic curves without singular points, such as the curves B for linv^oo Tfc jT .. 
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FIG. 6. As in Fig. H for k = (a) 5 (b) 6. Chromatic zeros for r = (a) 10 (b) 8 are shown for comparison. 



VIII. THE FAMILY S k>r = HEG k - 2 {K 3 + T r ) 



To illustrate one infinite family of graphs with noncompact B of the type HEG^i{K v + G r ) with p > 
we discuss an infinite family constructed by homeomorphic expansion of the bonds of the subgraph G r = 
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in the graph (K3 + T r ). The beginning family, K3 + G r , is of type (3.3) (and can be obtained by complete 
removal of all three of the bonds in the if 3 subgraph in K% + G r ). We denote this family as 



The number of vertices is given by 



Sk,r — HEGk-2\Ks + T r ) 



v(S ktr ) = (r-l)(k-l) + 4 



(8.1) 



•2) 



An example of a graph of this type is shown in Fig. ^|(c). The chromatic number and girth are the same 
as those of the graph T k ^ r , given in eqs. ( |6.6| ) and ( |6.7| ). By the deletion-contraction theorem, we find the 
recursion relation 



P(S k , r , q)= (q- 3)D k + (-If" 1 P(5 fc , r _i, q) + 3D fe P(T fc , r _ 1; q) - q{q - l)D k (D k+1 ) r 



(8.3) 



Solving this, we obtain the chromatic polynomial 



P(£fc,r,3) = 3 



q(q-l)D k+1 -P(T ka ,q) 



ar 2 - a r - 2 



+ q{q- i)ai 



or 2 - or 2 



+a r 3 - 2 P(S k , 2 ,q) 



where a\ and a 2 were given in eqs. (3-11) and (3.12), 

a 3 = (q-3)D k + (-l) k - 1 , 

and 

P(S ka ,q) = (q- 2)P(T ki2 , q) + q(q - l) 3 D k ^ 



(8.4) 



(8.5) 



(8.6) 



In the limit r — > 00 with k fixed |3^] , the locus B is determined by the degeneracy of magnitudes of leading 
terms 

M - |a 3 | (8.7) 

(02 is never a leading term in this case). The degeneracy equation in the y variable takes the form 

|1 - 2y + 3(-l) fe +V I = |1 + (-l) k+1 y k \ (8.8) 

Since y — is a solution of this equation, B is noncompact in the q plane, passing through y = z — 0. In 
polar coordinates, with y — pe l/3 , eq. (|Q|) yields 



2p 2k - 1 +p-cos/3+ (-l) fe -V i cos(kf3) - 3 /9 cos((yfc - 



= 



(8.9) 



As p — > 0, it follows that cos (3 = 0, i.e. /3 = ±7r/2, so that S approaches y = z = vertically 

To calculate the point at which the boundary B crosses the real q, or equivalently, y or z axes, we set 



/3 = in eq. (8.9); for p ^ this gives 

2p 2fe - 1 + p - 1 + (-l) fc -V fc_1 (l - 3p) = 



(8.10) 
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For both even and odd k, eq. ( 8.1 0| ) has only one acceptable (real positive) root for p, which is thus y c 



For odd k this root is y c = 1 (equivalently, z c = 1/2, q c = 2). For even k, the real positive crossing point 
increases monotonically from y c — 0.647799 (z c — 0.39313) for k — 4 to y c = 1 (z c = 1/2) as k goes to 
infinity through even values; i.e., q c decreases monotonically from 2.54369 to 2 over the same range. 

In Fig. [?] we show the respective boundaries B in the z plane for k — 3,4 |33|. Note that for k = 3, B 
is close to being, but is not, a circle. Again, for comparison, we show chromatic zeros for reasonably large 
finite values of r in each of these figures. As is evident from Fig. ^ and from the figures for higher values 



of k 33 1, B for lim r ^oo Sk, r displays some similarities with that for lim r ^oo Tfc jr . Hence we do not show 
these figures for higher values of k here. We remark, however, that there are some differences; for example, 
the number of connected components, N comp , is equal to the respective values 1,1,3,3,3,5 for 3 < k < 8 for 
linv^oo Tk, r but is equal to 1,1,3,3,5,5 for the present case, lim r ._ >00 Sk >r with 3 < k < 8. As was true with 
N C omp. for hnij.^oo X^r, for k > 4, the values of N comp , are only weakly bounded by the Harnack theorem. 

For each value of k, the region R\ is the one occupying the non-negative interval of the z axis adjacent to 
the origin. The region i?2 occupies the rest of the z plane for 2 < k < 4, and the rest of the z plane except 
for the interiors of the other closed loops for k > 5. We find 

W( lim S Kr , q) = (ax) 1 '**-!) = ODfc+r) 1 ^- 1 ) for q G R, (8.11) 



and 



\W([hm S kir iq)\ = \a 3 \ 1 ^ k - 1 ^\(q-3)D k + (-l) k ~ 1 \ 1 ^ k -^ for q G R 2 (8.12) 



For k = 5, 6 (k — 7, 8) there is one pair (are two pairs) of complex-conjugate closed loops disconnected from 
i?i , which we denote Rdiscj and R* disc j ■ In these regions we find 

\W([]im S k , r ],q)\ = (ai) 1 ^" 1 ) = {Dk^) 1 '^ for q e RdUcj , R* disc < (8.13) 

The feature observed in the r — > oo limits of the families T k , r and C k , r that the disconnected regions Rdisc,j 
and R* disc j tend to lie along the line z = 1/2 is also present in this case. 
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FIG. 7. Boundary B in the z = l/q plane for the r — > oo limit of the family of graphs HEGk-2(Kz + T r ) with k = (a) 3 
(b) 4. Chromatic zeros for r = 14 are shown for comparison. 



IX. GENERAL CONDITION FOR NONCOMPACT B 



The key ingredient to construct families of graphs with noncompact W boundaries B in the q plane and 
resultant reduced functions W re d. ({G}, q) that are nonanalytic at z = l/q = is to produce a chromatic 
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polynomial with the feature that there are two leading terms with a degeneracy condition (2.6) that has a 
solution at z = (cf. eq. (|2.9| )). A general statement of this condition was given as the theorem of Section IV 
of Ref. |l5| ]. From our studies in the present paper, we can add some geometrical comments to this algebraic 
condition. 

An important property of our families of graphs with noncompact W boundaries is that none of these 
families is a regular lattice graph. This is in accord with the derivation of the large-g expansion for regular 
lattices ^5j,|(|. (Of course, the property that {G} is a regular lattice is not a necessary condition for the 
associated boundary B to be compact in the q plane and hence W r ed.{{G},q) to be analytic at z = 0.) A 
basic feature of a regular lattice graph is that (except for boundary vertices, which yield a negligible effect 
in the thermodynamic limit and are absent if one uses periodic boundary conditions) all vertices have the 
same coordination number (— degree A in usual graph theory terminology). A generalization of this is 
encountered in duals of Archimedean lattices, such as the [4 • 8 2 ] (union-jack) lattice, where the vertices fall 
into a finite number of sets with different coordination numbers (e.g. A = 4 and A = 8 for the [4-8 2 ] lattice). 
In both of these classes of lattices, there are no vertices with infinite degree. 

In contrast, a property of the families of graphs that we have studied with noncompact loci B(q) in the 
limit r — * oo is that in this limit they all contain an infinite number of different, non-overlapping (and 
non-self-intersecting) circuits, each of which passes through at least two fixed, nonadjacent vertices. This 
immediately implies that these aforementioned nonadjacent vertices have degrees A that go to infinity in 
this limit. In the families that we constructed earlier fll5f and some of the additional families discussed in 
sections II and III of the present work, this property was produced by starting with a family of the form 
K p + G r and removing one or more bonds from the K p subgraph, thereby rendering two or more vertices 
that were originally adjacent no longer adjacent. As we showed in Section IV, it can also be produced by 
homcomorphically adding one or more vertices to a bond in the K p subgraph, thereby again rendering two 
vertices that were originally adjacent no longer so. It should be mentioned that the condition of having 
two or more vertices with linv^oo A(r) = oo is not, by itself, sufficient to produce a noncompact B{q). For 
example, consider the family of 11 p- wheels" , 

{Wh)W=K p + C r (9.1) 

(where n = p + r is the number of vertices) that we constructed in Ref. O] . The degree of each of the p 
vertices in the K p subgraph of the graph K p + C r is A = p — 1 + r, so that in the limit r — > oo (with p fixed), 
this degree A — * oo. However, the corresponding boundary B is compact; specifically, we showed that it is 
the unit circle \q — (p+ 1)| = 1 fl4f| . Furthermore, we observe that the presence of non-adjacent vertices with 
degrees A such that linv^oo A = oo is also not, by itself sufficient to guarantee that B(q) is noncompact. 
This is easily seen by considering, for example, two p = 1 p-wheel graphs, (W7i)l whose central vertices 
are connected by a tree graph containing at least two bonds. As n — * oo, the degrees A of the central 
vertices go to infinity, but again the resultant W boundary B is compact; indeed, it is the same as for a 
single p = 1 p-wheel, viz., \q — 2| = 1. Thus, a necessary feature is observed to be the presence, in the limit 
as the number of vertices of the graph goes to infinity, of an infinite number of different, non-overlapping 
(non-self-intersecting) circuits, each of which pass through the two or more nonadjacent vertices. 
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X. CONCLUSIONS 



In this paper we have further explored a fundamental problem in statistical mechanics - nonzero ground 
state entropy using as a theoretical laboratory the (/-state Potts antiferromagnet. We have presented 
a number of exact calculations of the zero-temperature partition function Z(G, q, T = 0) or equivalently 
the chromatic polynomial P(G,q), and the corresponding limiting function representing the ground-state 
degeneracy, ^({G}, q), for this model on various families of graphs G for which the boundary B of regions 
of analyticity of W in the complex q plane is noncompact, passing through z = 1/q = 0. The study of 
these graphs thus gives insight into the conditions for the validity of the large-g Taylor series expansions of 
the reduced function W red , ({G}, q). In addition to families obtained by removal of bonds from nonadjacent 
vertices in the K p subgraph of K p + G r , we have constructed and investigated a number of families of graphs 
by the powerful method of homeomorphic expansion from respective starting families. We have shown how 
the families thus obtained have, in the limit of an infinite number of vertices, noncompact boundaries in the 
q plane which pass vertically through the origin of the z — 1/q plane and have support for Re{q) > 0. 

This research was supported in part by the NSF grant PHY-97-22101. 



XI. APPENDIX 1 



In this Appendix we gather together some convenient formulas concerning the function Dk, defined in eq. 



( 1.19 ) in the text. One has 

D k {q = Q) = {-l) k {k-l) (11.1) 

and 

D k (q = l) = (-l) k (11.2) 

Since for k even, the circuit graph C k is bipartite, which is equivalent to the fact that the chromatic number 
x(Ck even) = 2, and thus P(C k even,q = 2) = 2, it follows that 

D keven {q = 2) = 1 (11.3) 

Since for k odd, x(G fe odd) = 3 and P{C k odd, q = 2) = 0, we have 

D k odd (q = 2) = (11.4) 

This zero results from a linear factor, i.e., 

D kodd (q = 2) = (q-2)Pol(q) (11.5) 

where Pol(q) is a polynomial of degree k — 3 in q with Pol(q = 2) =/= 0. Two identities that we have derived 
and used for our calculations are listed below: 
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Dk-aDk-! = (-l) fe (11.6) 

D k+1 - D k = D 3 D k + (-l) fc+1 (11.7) 
The proofs follow immediately from the definition of Dk- 
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